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INTRODUCTION 


This  is  the  last  of  several1-4  reports  which  conceptually  go 
beyond  the  usual  practical  treatments  of  Kfc  in  the  theoretical 
•  applied  mechanics  literature.  That  literature  consists  of 

derivations  of  a  Kk  or  other  single  associated  measure  (COD) ,  for 

» 

each  of  many  configurations  of  cracked  material  and  loading  such 
that  evaluation  of  the  appropriate  formula  to  any  test  under  Kk 
conditions  gives  the  same  toughness  value  as  would  a  similar  test 
of  any  other  configuration  of  the  same  material.  This  is 
consistent  with  the  notion  that  toughness  is  a  unique  fundamental 
property  of  materials,  one  not  made  up  of  simple  fundamental 
quantities  such  as  fracture  stress  or  work  hardening  modulus. 
Here,  however,  Kk  toughness  is  not  considered  basic  but  dependent 
on  two  metallurgically  controllable  parameters,  the  tensile 
brittle  fracture  stress  and  a  postulated  strain  parameter,  the 
effective  crack  tip  radius  existing  at  the  Kk  fracture  load 
stress;  this,  in  turn,  is  considered  to  be  dependent  on 
metallurgy,  e.g.,  sulfur  content,  as  well  as  the  Klc  loading 
stress.  The  reports  describe  how  to  determine  the  fracture 
»  stress  and  the  radius  from  Kk  tests  using  a  postulated  fracture 

criterion  characteristic  of  the  material  and  a  curve  derived  from 
the  solution  for  the  state  of  stress  in  the  yielded  region  of  a 
short  elliptically  shaped  flat  crack  lying  across  the  center  of  a 
plate  under  uniform  tensile  loading.  Here  a  more  exact  solution 
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supplants  the  approximate  ones  used  previously  and  it  is  shown 
that  radius  determinations  remain  unchanged  although  fracture 
stress  determinations  are  75%  as  large.  The  state  of  application 
and  advantages  of  the  fracture  stress  approach  over  the  energy 
approach  are  discussed. 

The  title  of  this  report  was  first  used  on  a  report4  issued 
in  October  1978;  for  comparison  that  report  shall  be  called  the 
First  Report.  It  was  a  "how  to"  report  based  upon  a  curve 
obtained  by  an  approximate  analysis  employing  elasticity  shear 
stress  trajectories  with  the  hope  that  the  analysis  would  soon  be 
followed*  by  a  more  exact  finite  element  one  s-6.  The  latter  was 
expected  to  apply,  not  only  to  the  crack  problem,  but  to  evaluate 
the  approximate  type  of  analysis  for  the  crack  and  other  stress 
concentration  problems  where  the  yield  point  was  locally  exceeded 
and  for  which  an  exact  analysis  was  difficult  to  achieve. 

A  careful,  time  consuming  finite  element  analysis  was 
finally  made  in  1984  by  Colin  E.  Freese  and  Dennis  M.  Tracey 
largely  in  answer  to  the  above  expectation7'10.  It  is  the  one  used 
in  this  report,  is  considered  very  good  (the  best  available) ,  and 
contains  a  surprising  limiting  stress  result.  However,  it  leaves 
some  important  questions  unanswered:  whether  the  stress  is 
maximum  just  off  the  center  line  of  the  crack,  as  found  by  the 


it  vu  followed  by  *  mandating  finite  element  and  experimental  aoatyeif  led  by  A  fecit  Kabtyaafai  (Reference*  3,6).  The  experimental  rendu  were 
very  util  factory,  but  it  was  evident  feat  a  finer  grid  would  be  neceaaary  for  preciae  finite  dement  reaulta. 
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approximate  method  (a  distance  too  small  for  the  finite  elements 
used)  and  whether  it  is  appropriate,  even  near  the  crack  tip,  to 
compare  the  approximation  which  was  based  upon  the  uniform 
biaxial  loading  used  by  fracture  mechanics  pioneers,  following 
Griffith11, ",  with  the  uniaxial  tension  loading  used  by  Freese  and 
Tracey.  The  delay  in  publishing  the  present  analysis  was  largely 
due  to  the  expectation  that  time  could  have  been  found  for 
definitive  answers  to  such  questions. 

The  approach  here  to  crack  toughness  or  brittleness  is  not 
an  energy  approach  but  a  tensile  brittle  fracture  stress14 
approach  (Fig.l).  Using  shear  stress  trajectories  (Fig.  2),  the 
maximum  tensile  stress,  Sa  for  loading  stress  S,  is  computed  from 
Sm  =  Y  (1  +  A)  where  Y  is  the  plane  strain  yield  strength  and  L 
is  the  largest  angle  change  found  among  all  the  principal  shear 
stress  trajectories  crossing  the  yield  region  at  the  crack  tip. 
Thus,  this  stress  is  ahead  of  the  crack  and  increases  with  the 
amount  of  yielding;  hence,  load;  if  it  equals  the  brittle 
fracture  stress,  F,  of  the  material,  failure  occurs.  The  strain 
where  the  trajectory  corresponding  to  L  crosses  the  yield  region 
is  small;  hence,  work  hardening  is  negligible,  and  since  the 
trajectory  is  far  from  the  crack  tip  with  any  appreciable  L,  the 
exact  shape  of  the  tip  is  of  limited  importance;  the  radius  of 
the  tip  found  from  the  analysis,  whatever  its  origin,  is  an 
effective  radius. 
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The  finite  element  solution  is  based  upon  a  crack  tip  radius 
of  p  =  0.001,  taking  half  the  crack  length  as  the  unit  of 
distance,  and  was  applied  to  other  cases  through  dimensionless 
ratios.  For  that  solution,  the  dimensionless  ratio  L  was 
computed  by  equating  Y  (1  +  L)  to  the  maximum  stress 
corresponding  to  an  applied  loading  stress.  Y,  the  plane  strain 

yield  strength,  was  taken  to  be  2xl05/v'3. 

FINITE  ELEMENT  DATA  AND  MATCHING  FORMULAE 

Our  notation  and  definitions  are  contained  in  Table  I. 

The  finite  element  program  models  a  long  flat  tensile 
specimen  which  is  loaded  uniformly  at  the  ends  and  is  pierced  by 
a  central  crack  of  length  2a  across  its  middle. 

For  present  use,  Freese  supplied  the  unsmoothed  finite 
element  method  results,  as  well  as  L,  M,  and  L  directly  computed 
from  it,  by  Table  II.  The  L  and  M  versus  L  data  is  used  to 
determine  fracture  stress  F,  and  effective  crack  tip  radius,  p, 
from  appropriate  (pressure  free  crack  surface)  experimental  K,c,  Y 
(toughness,  yield  strength)  data. 
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TABLE  I 

NOMENCLATURE  AND  BASIC  RELATIONS 

a  =  crack  length/ 2 

p  -  0.001  *  crack  tip  radius  in  Table  II 

Yt  =  Tensile  test  yield  strength  =  102  ksi  in  Table  II 

Y  =  2YcA/3  =  Plane  strain  yield  strength,  the  yield 
strength  used  in  all  analyses  in  this  report. 

S  =  The  uniform  tensile  loading  stress 

x  =  Distance  from  crack  tip  to  point  of  maximum  stress,  Sm 

Sm  =  Maximum  stress  ahead  of  crack  at  loading  stress  S, 
i.e.,  stress  at  x. 

F  =  Tensile  stress  needed  to  cause  fracture;  here  the  nil- 
ductility  tensile  fracture  stress. 

Sm  =  F  under  K,c  test  conditions. 

L  =  Maximum  shear  stress  trajectory  angle  change; 

Sm  =  Y(1  +  L) 

L  =  Y/Syfap 

=  Yjnp/K,  a  dimensionless  parameter 
=  3.65148/S  in  Table  II;  at  initial  yielding  L  =  2 

K  =  Sy/na  here,  =  KIc  toughness  under  K,c  test  conditions 
M  =  1a[xT P  -  {Y/K)yJ%x,  by  definition. 
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TABLE  II 

DIMENSIONLESS  QUANTITIES  COMPUTED  FROM  FREESE-TRACEY  DATA  FOR 
AN  ELLIPTICAL  CRACK  IN  A  TENSILE  SPECIMEN 


a  j 


MATCHING  FORMULAE 


The  following  L  and  M  formulae  are  selected  from  the  many 
that  were  devised  to  match  and  smooth  the  L*  and  M  data  tabulated 
in  Table  II  and  shown  in  Figures  3  through  5. 


With  F  =  Y(l+A.)  for  no  pressure  on  crack  faces, 

u.iTSp/K)  =  T35- 1"  TTirk  ;  0si41-11 


1 

1.29 


In 


1.46 

SjY- 0.88 


when  F=  F  ( 1  +  A.) 


M  (  » UJx/ p  a  Yjnx/K)  =  Ltsinh  [  (8/3)  A.]  /2)'/t 


=  L(sinil  [  (8/3)  ( Sm/ Y- 1 )  ]  / 2},/4 


L  maximum  is  2,  and  goes  precipitately  to  effectively  zero  at 
A.-1.11,  0  sA.sl.ll 


x/p  maximum  is  5,  0sx/p<;5 


Sm/Y  maximum  is  *  2.11Y 


Considering  the  validity  of  dais.  etc.,  we  actually  use  the  formula  for  L  to  d  -  1.11  =  10/9,  where,  in  fact  L  seems  to  sharply  drop  to  zero; 
approximately  tee  this  drop,  multiply  the  L  formula  by  R=|l-(8d/9)"p ".  The  Tables  for  L  vs.  d  of  the  First  Report  may  be  used  in  place  of  the  above 
vs.  d  formula  if  the  d’s  are  multiplied  by  3/4;  see  L,  vs  d,  p.  17. 
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If  there  is  plane  strain  pressure,  P,  on  the  crack  faces13, 

sa  =  y(i+4.)  -p. 

The  data  is  not  precise  below  L  =  0.25,  say,  a  caution, 
especially  for  those  who  use  the  data  directly  rather  than 
through  the  formulae. 

BRITTLE  (MIL-DUCTILITY)  TENSILE  FRACTURE  STRESS , 
EFFECTIVE  RADIU8,  YIELD  8TREM6TH 

FRACTURE  STRESS 

If  the  tensile  test  temperature  of  a  structural  metal  is 
lowered  sufficiently  the  metal  is  found  to  break  with  little  or 
no  ductility.  This  stress  is  our  fracture  stress,  F,  if  there 
has  been  slight  straining,  but  not  necessarily  if  the  temperature 
is  so  low  that  there  is  no  strain  for  there  is  evidence  in  at 
least  some  cases,  that  a  completely  undeformed  specimen  requires 
a  higher  stress  for  fracture  than  a  slightly  deformed  one.  In 
our  case  the  greatest  stress  lies  on  the  border  of  deformed 
material  beyond  the  front  of  the  crack.  See  Figure  6,  where, 
assuming  uniform  yielding,  F  is  about  60  ksi  and  high  strain  rate 
takes  the  place  of  low  temperature. 

F  is  of  course,  a  macroscopic  brittle,  essentially  nil- 
ductility  tensile  fracture  stress,  not  the  theoretical  fracture 
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stress  of  a  homogeneous  flawless  material;  the  theoretical  stress 
is  presumably  reduced  to  it  by  heterogeneity,  flaws  or 
inclusions. 

For  our  conception  of  fracture  stress  see  Fig.  7. 

In  the  next  section  on  use  of  the  finite  element  data,  we 
assume  that  for  certain  ranges  of  data  that  a  linear  fracture 
stress  law  may  be  used,  i.e.,  Y  =  (1  +  0)Y  +  f,  where  0  and  f  are 
constants . 

For  the  high  strength  steels  with  which  he  was  familiar,  the 
author  observed,  or  inferred,  that  for  tests  at  low  to  moderate 
temperatures  covering  practical  ranges  of  brittle  failure,  that 
when  the  temperature  of  the  tests  alone  was  varied  F  =  f,  a 
constant*,  and  that  when  the  tempering  temperature  of  the  test 
specimens  alone  was  varied  F  =  Y  +  f,  i.e.,  the  fracture  stress 
was  a  constant  amount  above  the  yield  strength,  Y.  These 
relations  he  generalized  into  a  linear  fracture  stress  "law" 
which  he  assumed  suitable  for  the  range  of  stress  involved  in  the 
K,e  determinations  of  p  and  F,  i.e., 

F  —  (1+0)  Y  +  f 


Accordingly  the  writer  observed  dial  failure  in  guru  occurred  when  the  yield  strength  was  raised  by  low  temperature  operation  so  that  the 
concentrated  yield  strength  equalled  P,  the  fracture  strength  of  the  material  of  which  they  were  made.  (See  Discussion) 


9 


Thus  /3  =  -1  for  the  varied  test  temperature  tests  above  and 
-  0  for  the  varied  tempering  temperature  tests,  above. 

RADIUS 

In  the  usual  small  deformation  elastic  solution 
corresponding  to  a  flat  elliptical  notch  in  a  tension  field,  the 
tip  radius  increases,  in  proportion  to  the  load,  depending  on  the 
modulus  of  elasticity,  E;  this  is  likewise  the  case  for  the  non- 
cyclic  loading  elastic-plastic  solution  except  that  the 
proportion  depends  on  the  amount  of  yielding  and  thus  on  the 
yield  strength  and,  to  a  more  limited  extent,  on  the  work 
hardening  modulus.  The  point  here  is  that  there  is  a  unique 
radius  for  each  load;  the  metallurgically  dependent  radii  we  have 
postulated  do  not  exist  for  that  fixed  geometry. 

But  our  deductions  of  radii,  based  on  the  curve  fitting 
methods  described  in  the  next  section  are  evidence  that  a  change 
does  take  place  at  the  very  highly  stressed  tip,  that  we  get 
radii  other  than  above. 

And,  Taggart'3,  fitting  a  parabola  to  the  deflection  of  crack 
boundaries  near  the  tip,  found  an  increase  in  p  with  K,  an 
increase  to  KIc  load.  The  increase  was  not  the  same  for  different 
materials,  as  it  would  be  if  the  crack  tip  radius  increased 
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uniquely  with  load  for  a  given  Y  (and  work  hardening  modulus,  if 
large)  . 

Thus,  although  F  varies  with  Y,  it  is  assumed  in  the 
analysis  for  p  below,  that  p  remains  the  same  in  any  two  ordinary 
test  conditions  for  Kfc  unless  the  fracture  shows  a  change  such 
as,  e.g.,  the  appearance  of  cleavage  fracture  in  the  one  test  but 
not  the  other. 

Theoretically,  if  not  practically,  according  to  the 
definitions  of  L  and  M,  x  the  distance  to  the  point  of  maximum 
stress  along  the  crack  axis  could  have  been  chosen  instead  of  p 
as  one  of  our  crack  characterization  parameters. 

YIELD  STRENGTH 

Perhaps  this  is  as  good  a  place  as  any  to  point  out  that  a 
likely  source  of  error  is  the  use  of  conventional  determinations 
of  yield  strength  in  our  formulae.  Our  yield  strength  should  be 
that  of  completely  yielded  material  such  as  that  determined  by 
cyclic  tests,  even  if  the  yielding  is  small.  The  conventional 
determinations  are  generally  too  high16'18;  a  particularly  bad  case 
is  that  of  mild  steel  in  which  even  the  lower  yield  point 
represents  non-uniform  yielding,  where  the  first  yield  is  that  of 
a  slip  band  at  an  enormous  rate  of  strain  [(specimen  length/slip 
band  width)  x  nominal  strain  rate].  It  is  arguable  that  the 
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endurance  limit  stress  is  the  proper  one;  taking  the  yield 
strength  to  be  half  the  tensile  strength  may  be  acceptable  in 
some  cases. 

The  fracture  stress  may  be  little  above  the  apparent  yield 
strength,  but  substantially  above  the  actual  yield  strength. 

USE  OF  THE  FINITE  ELEMENT  DATA 

MODE  I  STANDARD  TEST  DATA 

The  object  here  is  to  determine  F  and  p  from  Kfc  data 
obtained  with  pressure  free  crack  faces.  If  only  L  (=  Yjnp/Klc) 

is  known  for  two  or  more  cases,  L  must  be  supplemented  by  other 
information  such  as  an  appropriate  fracture  law,  here  taken  to  be 
F=(l  +  j8)Y  +  f,  and  assumptions  supported  by  metallurgy  and 
fractography  of  constancy  of  p  and  F  in  the  range  of  data  used  in 
the  determination. 

Examples  will  be  found  in  the  First  Report  from  which  we 
have  drawn  Fig.  8.  From  this  we  note  how  well  F  and  p 
determinations  correspond  to  the  experimental  data  as  well  as  the 
effect  of  Sulfur  level  on  the  radius.  Four  cases  are  cited  here, 
defined  by  what  is  used  and  constant  (C  *  constant)  in  the  range 
considered. 
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M  with  L 


F  and  p  can  be  determined  from  the  equations  for  L  and  M  if 
x  can  be  measured  or  inferred,  e.g.,  from  ripples  which  might 
result  from  successive  fracturing  and  opening  to  the  crack  tip, 
especially  if  the  incipient  fracture  is  just  off  the  center  line 
of  the  crack,  x  being  the  distance  from  the  tip  of  the  crack  to 
the  site  of  incipient  fracture,  x  may  be  of  the  order  of  one  to 
five  thousandths  of  an  inch. 

Thus  in  M  we  would  have  an  expression  in  which  all 
quantities  except  F  are  known;  thus  M  determines  F.  Insertion  of 
this  F  into  L  is  an  expression  for  p.  A  hand  held  computer  which 
is  programmable  is  handy  for  this. 

L,  0  =  -1,  C(p,F) 

For  0  =  -1,  corresponding  to  constant  F  at  different  test 
temperatures,  F  =  f ,  a  constant  though  Y  changes.  (See  Fig.  9) . 
Assuming  F  and  p  are  the  same  for  two  data  points,  Y,  K,e,  two 
simultaneous  equations  may  be  formed  and  solved  for  F  and  p  since 
each  insertion  of  a  data  pornt  Y,  KIc  into  the  expression  for  L 
provides  an  equation.  Obviously,  different  combinations  of  data 
points  may  be  used  to  check  a  result. 
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L,  0  a  0,  C(p,f) 


For  0  =  0,  corresponding  to  tests  at  different  tempering 
temperatures,  F  =  Y  +  f  where  f  is  a  constant  and  Y  changes.  (See 
Fig.  8)  Since  we  assume  p  -  constant  (subject  to  constancy  of 
existence  of  the  weakest  kind  of  fracture,  e.g.,  grain  boundary 
fracture) . 


L*Y^np/K  -  function  (F/Y) 

we  have  for  two  data  points: 

Yl}/np /Kx  -  functi ond+f/Yj 

and 

Y2yfnp/K2  -  function  (l+f/Yz) 

These  may  be  solved  simultaneously  to  find  f  and  p. 

L,  0,  c  (p,  F)  and  L,  C  <p,  f) 

The  scatter  in  experimentally  determined  values  of  K  is 
often  so  pronounced  that  a  curve  fitting  procedure  using 
transparent  overlay  graphs  and  involving  more  Kfc,  Y  values  is 
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recommended,  as  described  in  the  First  Report  and  repeated,  in 
part,  in  Appendix  I  -  Overlays.  (See  Figs.  10-12.) 

It  will  be  recognized  by  some  that  this  overlay  procedure 
can  be  done  very  readily  on  a  computer;  for  example,  following 
the  procedure  of  Kfe  versus  Y  for  various  tempering  temperatures 
of  steel  for  which  0  =  0,  plot  In  (1/AX)  versus  In  (1/A.)  from  the 
equation  for  L  versus  A.,  the  overlay,  and  separately,  plot  the  Klc 
versus  Y  data  as  In  K  versus  In  Y.  Move  the  latter  plot  of 
points,  along  with  its  coordinates,  parallel  to  the  coordinates 
of  the  In  (1/AX)  versus  In  (1/A.)  plot  until  the  points  fit  the 
curve  of  the  latter  plot.  Then  f,  a  constant,  is  that  coordinate 
value  of  Y  where  the  In  Y  coordinate  is  crossed  by  the  unit 

vertical  coordinate  of  the  overlay  and,  similarly,  f^/np  ,  a 

constant,  is  that  coordinate  value  of  K  where  the  In  K  coordinate 
is  crossed  by  the  unit  horizontal  coordinate  of  the  overlay. 

Then  F  =  Y  +  f  where  Y  is  any  Y  and*  p  is  {f^/np)2/nf2. 

Generally,  a  master  reusable  In  In  plot  of  [l/(A.-jS)L]  versus 
[1/(A.-0)J  is  made  on  transparent  In  In  paper  as  a  bank  of  /3 
constant  curves;  this  overlay  is  placed  over  any  particular  plot 
of  K,  Y  data  on  the  same  size  of  In  In  paper  as  the  transparency 
and  slid  parallel  to  the  coordinates  to  the  &  curve  which  best 
matches  the  data;  this  determines  0  and,  simultaneously,  f  and  p; 


From  experience  it  it  eeUom  worthwhile  to  allow  for  the  eltatic  chance  in  p  aa  ia  done  in  the  Fint  Report. 
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f  is  that  coordinate  value  of  Y  where  the  Y  coordinate  is  crossed 


by  the  unit  vertical  coordinate  of  the  overlay  and  fjn p  is  that 

coordinate  value  of  K  where  the  K  coordinate  is  crossed  by  the 
unit  horizontal  coordinate  of  the  overlay. 

HYDROSTATIC  PRESSURE  MODIFICATION  OF 

Kfc  formulae  obtained  by  energy  methods  are  independent  of 
hydrostatic  pressure1.  This  particular  consequence  may  be 
justified  by  the  belief  that  there  is  always  an  infinite  stress 
singularity  at  the  tip  of  the  crack. 

Here,  however,  it  is  believed  that  for  metallic  materials  of 
construction,  at  least,  there  is  an  effective  crack  tip  radius  at 
any  load;  that  consequently  all  stresses  are  finite  and  reduced 
by  the  addition  of  hydrostatic  pressure;  that  Kfc  failure  will 
occur  if  the  largest  stress  equals  the  tensile  brittle  fracture 
stress,  F,  of  the  material. 

Thus  we  now  consider  how  to  modify  Kfc  determinations  to  make 
them  applicable  to  the  case  where  there  is  uniform  pressure  on 
the  crack  surface,  a  pressure  that  might  be  applied  by 
superimposing  on  the  standard  loading  a  hydrostatic  pressure. 
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A  plot  of  L  vs.  L  for  a  crack  in  a  distant  field  of  uniform 
tension  S*  is  a  plot  of  l( S/Y)  (a/p)14]'1  vs.  L  where  if  S/Y  and 
a/p  are  known  one  finds  L.  The  addition  of  a  plane  hydrostatic* 
pressure  P  does  not  change  the  yielded  areas,  their  boundaries  or 
the  geometry;  it  only  makes  the  size  slightly  and  negligibly 
smaller.  It  does  not  affect  the  yield  boundaries  of  the  lobes  of 
yielded  regions  at  the  tips  of  the  crack;  therefore  L  is  not 
changed.  P,  however,  is  added  (negatively)  to  every  component  of 
normal  stress,  including  boundary  loading  stress,  S,  but  we  still 
think  of  the  boundary  loading  stress  (S  -  P)  as  made  up  of  its 
two  components  since  only  S  controls  L;  L  in  L  vs.  L  does  not 
contain  P;  it  is  still  [( S/Y)  (a/p)14]'1.  Though  L  is  unchanged  it 
is  used  to  determine  S.;  if  there  is  no  P,  we  have  Sm  =  Y(1  +  A.) 
but  if  there  is  an  addition  of  hydrostatic  pressure  P,  as 
described***  above,  Sm  =  Y(1  +  A.)-P. 

Thus  without  P,  A.  =  Sm/Y  -  1  and  with  P,  A.  =  (Sm  +  P) /Y  -  1 
although  K  is  independent  of  P  in  L  =  Yy/np/Sy/na  =  Yy/np/K. 


The  Griffith  cue.  but  equivalent  For  an  ideal  crack  to  the  cue  of  uniaxial  tension  perpendicular  to  the  crack. 

•• 

Obviously,  the  orifinal  loadinf  could  have  been  in  the  crack  inalead  of  applied  externally  in  which  case  P  (if  P  was  the  loading)  would  be  in  L 
instead  of  S.  We  choose  the  external  case  here  because  it  is  the  case  for  which  K,,  measurements  see  made. 

••• 

According  to  the  equilibrium  equations,  for  a  material  where  it  undergoes  a  constant  yielding  stress  ,  Y,  A.  corresponds  to  so  increase  A.  Y 
over  the  magnitude  of  the  principal  stresses  at  the  crack  tide  of  the  boundary  of  the  yielded  region  (lobe  boundary)  by  the  corresponding  stresses  located  by 
following  a  principal  shear  stress  trajectory,  on  the  other  side.  Thus,  die  magnitude  of  a  principal  stress  (Y  -  P)  on  the  crack  side  becomes  (Y  -  P)  +  A.Y  = 
Y(1  +  A.)-P  for  a  corresponding  one  on  the  other  tide  u  the  difference  in  directions  is  A.  This  is  the  same,  of  course,  as  increasing  the  original  S_  stress, 
Y(l  +  A.)  by  the  hydrostatic  stress  -P. 
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Obviously,  if  S,,  =  F,  a  material  constant,  A.  is  greater  on 
the  L  vs.  A.  plot  with  P  than  without  P. 

Thus,  a  knowledge  of  F  and  p  found  from  a  standard  test 
specimen  Kfe  can  be  used  to  determine  what  'Kfc'  value  should  be 
used  if  there  is  known  internal  pressure  P:  A.  *  (F  +  P)/Y  -l, 
this  determines  L;  then  =  Yy/np/L,  calling  the  critical 
%Kfc'  with  pressure,  a  quantity  larger  than  Kfc. 

With  guns,  |P|  -  S,  so  there  is  no  external  loading  stress, 
only  an  internal  pressure  (S) . 

QUANTITATIVE  RELATIONSHIP  BETWEEN  THE  FINITE  ELEMENT 
AND  FIRST  REPORT  RESULTS 


L  versus  A. 

A  comparison  of  the  data  for  the  basic  L  versus  A.  curve  of 
the  First  Report  and  the  present  Freese-Tracey  finite  element 
data  indicates  to  possibly  95%  accuracy  that  for  any  specified 
value  of  L,  and  thus  any  combination  of  Y,  K,c  and  p  making  up 
that  value  of  L,  in  either  source, 

L  =  (3/4)  A., 

where  the  subscript  "1"  stands  for  the  First  Report.* 


That  the  L  vs  4.  Tables  of  the  First  Report  may  be  used  to  L.  cut-off  if  the  tabulated  6-’i  are  multiplied  by  3/4. 
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Since,  also 


F  *  Y(1  +  L) 

F,  -  Yt(l  +  Lx) 

we  have 

F/Y  -1  *  (3/4)  (F,/Y,  -1) 

where  L  =  Yy/np/KIC  =  Yiy/npl/ KIcl 

In  particular,  if  Y  =  Y,,  so  that  it  is  also  necessary  that 
Vp/Kic  =  JVJKIC1. 

F  -  Y  -  (3/4)  (F1-Y) 

which  may  be  rewritten  for  different  insight  as 

F  =  Y  +  (3/4)  (Fx-Y) 

F  =  3/4  Fx  +  Y/ 4 

F  -  Fx  -  (Fx  -  Y)  / 4 
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r  =  (i  +0)Y  ♦  f  v*.  a. 


Next,  we  note  how  the  linear  fracture  stress  laws  compare, 
also  using  F  -  K(l  +  A.)  and  Fx  =  J^d+A^)  . 

Since  we  wish  the  consequence  of  a  particular  assigned  value 
of  0  in  either  case,  0  =  0,;  in  designating  0  we  are  thus 
defining  F  in  terms  of  Y  and  f: 

F  =  (1  +p)  Y+f 

Fx  =  (1+  p  )  Y1+f1 

We  note  the  result  of  elimination  of  0  from  these  as  well  as 
two  cases  of  special  interest,  0=0  and  0  =  -1. 

Bv  elimination  of  8 

£/Y-£jYx  =  (F/Y-l)  -  {FjYx-D 

whence,  since  F  =  Y(1  +  L)  and  Fx  =  l^d  +  A^)  , 


£/Y-£jYx  = 


and  thus  if  Y  =  Y,, 
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f~fx  =  Y(A.-A.X) 


where  we  can  substitute  for  the  angles  their  expressions  in  L. 

For  8  =  0.  with  Y,  =  Y: 

F  =  Y+f,  SO  F/Y-l  =  A.  =  f/Y 

Fx  =  Y+/1#  SO  FjYx- 1  =  A.x  =  fx/Y 

A.  =  (3/4)  A.x;  f/Y  =  (3/4)  ( fjY ) 

i.e. ,  f  *  (3/4)  fj 

We  noted  in  L  vs.  L  that  L  =  (3/4)  A.,  used  here,  was  for 
equality  of  L's  in  the  expressions  for  A.  vs.  L. 

Thus, 

P  =  Pi 

since  we  assume  Y  =  Y,  and  KIc  =  KIcl  in  the  expression  for  L. 

See  Figure  11,  the  /3  =  0  overlay  for  the  First  Report,  on 
which  new  axes  have  been  added  making  it  valid  for  the  new  finite 
element  data  A.  =  (3/4)  A.,;  these  axes  are  displaced  downwards  and 
to  the  left  by  In  (3/4)  from  the  old  ones;  this  follows  from  the 
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fact  that  on  In  In  paper  the  old  coordinates  are  1M,L  vs  1  / Lx 
while  the  new  ones  are  l/[(3/4)AnL]  vs.  1/0/4)^,  i.e., 

on  In  In  paper  we  have  In  (1  / L^L)  -  In  (3/4)  vs.  lnl/A^  -  In  (3/4)  . 

For  B  =  -l  the  strength  law  gives 

F  =  f 
F,  =  f, 

From  this  because  A.  =  (3/4)  A.,  we  have  for  Y  =  Y,, 

F/Y-l  =  (3/4)  (FjY- 1) 


i.e. 


9 


F  -  (3/4)  Ft  +  Y/ 4 

as  was  expressed  above  under  L  vs.  L.  although  the  strength  law 
is  designed  to  make  F  independent  of  change  in  Y  for  j8  =  -1.  The 
consequences  are  discussed  below. 

Using  the  First  Report  data,  if  one  found  F  and  p  for  0  =  -1 
by  precisely  matching  experimental  observations  of  Klc  and  Y,  one 
could  not  get  a  precise  match  for  F  and  p  using  the  new  data  for 
0  =  -1;  one  would  have  to  change  /3  to  get  a  match;  in  this  case, 

0  =  -1,  would  change  to  0  =  -3/4: 
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1 


1 


1 

^-p 


(3/4)  +  3/  4  ‘  3/4  (A^  +  l) 


where  the  second  expression  is,  of  course  L  +  3/4.  On  In  paper 
the  3/4  is  simply  a  translation  movement;  no  change  in  shape 
takes  place.  But  this  is  not  the  desired  match  for  /3  =  -l  for 
the  finite  element  solution  if  the  match  was  made  for  0  -  -1  for 
the  First  Report. 

Thus,  believing  in  the  utility  of  the  /3  =  -1,  i.e.,  F  =  f, 
independent  of  Y  result  for  tests  at  various  test  temperatures, 
we  conclude  that  there  was  not  an  exact  correlation  for  0  =  -1 
tests  in  the  First  Report,  though  a  useful  one. 

What  we  find  is  that  a  very  satisfactory  correlation,  a 
satisfactory  match  for  the  useful  range  of  data,  exists  between 
the  j8  =  -1  overlays  for  the  new  and  old  data,  that  is  the 
overlays  whose  axes  intercepts  on  the  In  In  plot  of  K,c  vs.  Y 

determine  F  and  p.  For  the  First  Report  the  overlay  is  a  In  In 

plot  of  1/  [  (A.i  +  1)  L]  vs.  1/(A.1  +  1)  ;  for  the  new  work  it  is 

1/  [  (0 .75^  +  1)  L\  vs.  1/ (0 .75^ +  1)  ,  since  A.  =  0.75^.  These  two 

plots  match  guite  well  if  the  horizontal  axis  of  the  latter  is 
displaced  to  In  0.935  downwards  and  the  vertical  axis  is 
displaced  leftwards  to  In  0.855,  of  the  original  First  Report 
axis. 
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These  new  axes  are  shown  in  Figure  12. 


The  intercepts  determining  F  and  p  are  displaced  the  same 
amount . 


Thus, 


F  =  0.855F! 


and 


F/np  -  0.935 


i.e. , 


0.855FlV/np  =  0.935Flvfifp7 


p  —  1.196  pj 


DISCUSSION  AND  CONCLUSIONS 


Before  taking  the  Freese  -  Tracey  data  as  definitive  and 
comparing  it  with  those  of  the  First  Report,  it  should  be  pointed 
out  that  the  maximum  stress  reported  in  the  First  Report  was  not 
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on  the  center2  line  of  the  crack  but  just  to  one  side  of  it 
(except  for  initial  yielding)  while  the  maximum  stress  reported 
by  Freese  -  Tracey  was  on  the  center  line;  the  grid  used  was  too 
coarse  to  make  such  a  distinction.  It  should  also  be  pointed  out 
that  Freese  -  Tracey  data  is  scattered  for  small  degrees  of 
yielding,  and  equal  stress  loading  along  and  perpendicular  to  the 
crack  was  used  in  the  First  Report  and  simple  tension  in  the 
Freese  -  Tracey  case  (a  difference  possibly  of  consequence  if 
yielding  is  relatively  extensive) .  The  shear  case  was  solved 
(not  published)  by  the  method  of  the  First  Report  and  showed  no 
concentration  of  stress;  a  result  which  seems  to  conform  with 
experiment,  thus,  the  method  is  of  interest  in  considering  the 
above  points. 

Except  for  a  sort  of  scaling  factor,  correlations  of  Klc,  Y 
data  to  find  F  and  p  by  Freese  -  Tracey  and  the  First  Report  are 
practically  equivalent.  Thus,  to  use  the  excellent  correlations 
in  the  examples  of  the  First  Report,  one  only  need  to  change  for 
0  =  0,  f  =  (3/4)  f|,  p  =  p,  and  for  0  =  -1,  F  =  .855  F,, 
p  =  1.196  p,  where  the  M,"  refers  to  the  First  Report  results. 

However,  there  is  a  very  important  difference  in  the  two 
sets  of  data:  there  is  an  upper  limit  to  L  found  by  the  finite 
element  approach,  L  =  1.125,  corresponding  to  a  limiting  distance 
of  yielding  ahead  of  the  crack  of  x/p  =  5.  This  angle  means  that 
the  greatest  stress  there  can  be  in  this  basic  crack  problem  is 
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(1  +  A.)  Y  *  2.11Y;  this  is  less  than  the  (1  +  n/2)Y  of  the  slip 
line  theory  for  the  straight  line  crack  and  much  less  than  might 
be  expected  from  the  First  Report. 

What  this  means  in  practice  is  that  if  the  fracture  stress 
of  the  material,  F,  is  greater  than  2.11Y,  no  Kk  crack  failure 
will  occur. 

This  kind  of  limitation,  a  characterization  which  seems  very 
important  in  practice,  is  even  stronger  than  was  anticipated  by 
former  analyses,  if  not  by  experiment  or  practice.  As  an 
example,  for  some  pieces  of  ordnance  equipment,  the  fracture 
stress  F  of  the  material  of  which  it  was  made  was  flat  over  the 
temperature  range  of  its  possible  employment;  when  the  equipment 
was  used  (i.e.,  loaded)  at  its  normal  temperature  of  operation 
there  was  no  brittle  Kfc  failure,  but  when  the  equipment  was 
loaded  at  a  reduced  temperature,  Y  became  sufficiently  great  for 
brittle  Kfc  failure  to  occur.  The  factor  of  concentration  was 
F/Y;  i.e.,  (l  +  A.)  See  Figure  9. 

Obviously,  the  fracture  stress*  of  the  material,  F,  should 
be  found  by  separate  experiment  if  possible.  Also,  as  far  as  is 
known,  the  two  physical  quantities  F  and  p  are  metallurgically 
independent.  If  this  is  so,  the  metallurgist  has  two  fabrication 

e 

Then  are  certain  metals  in  certain  stales  for  which  GrifTilh’s  surface  energy  density  seems  to  hold  using  his  formulation.  Il  is  doubled  chat  K  failure 
values  for  structural  metals  have  beta  appropriately  studied  for  comparison  or  otherwise  in  the  stale  in  which  they  are  macroscopicaily  brittle,  i.e.,  well  below 
the  transition  temperature,  a  state  usually  avoided  in  practice,  but  havnf  some  uses. 
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options  for  attaining  a  suitable  Kk.  There  is  evidence  in  the 
First  Report  that  p  is  dependent  upon  structure  and  sulfur  (see 
the  First  Report) .  If  F,  p  (possibly  machined  in)  and  Y  are 
specified  K  can  be  found  from  L  *  f  (A.) ,  i.e.# 

5V*p/*=f  (F/3T-1)  . 

In  the  case  of  the  example  above,  the  brittle  fracture 
stress  could  be  inferred  since  a  temperature  low  enough  for 
brittle  fracture  to  occur  in  a  simple  tension  test  could  be 
attained  using  available  coolants  such  as  liquid  nitrogen.  This 
is  the  case  for  some  other  materials  such  as  mild  steel  and 
tungsten  (see  Figure  6) ,  which  illustrates  stress-strain  curves 
with  their  fracture  envelopes.  Many  such  curves  were  drawn  to 
represent  the  results  of  many  tests  on  steel  at  this  laboratory. 
The  "nil  ductility"  brittle  fracture  stress  is  considered  to  be 
almost  nil  ductility  since  it  is  possible  that  a  test  under 
hydrostatic  tension  on  virgin  material  would  yield  a  higher  value 
corresponding  to  a  microstructure  unaltered  by  plastic  straining; 
this,  if  true,  is  of  no  concern  here  because  the  present  model 
presumes  yielding  and  fracture  at  the  boundary  of  plasticity. 

Other  attempts  to  determine  fracture  stress  have  been 
inconclusive  because  they  are  unfinished.  They  are: 
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•  Attempt  to  use  certain  notch  geometries  in  tension 
tests  with  partial  yielding. 

•  Attempt  to  put  uniform  exterior  tension  loading  on  a 
plate  containing  a  hole;  the  reason  for  using  this 
arrangement  is  that  with  it  there  is  a  maximum  stress 
which  occurs  at  the  boundary  of  the  plastic  region 
(which  extends  outward  from  the  bore)  and  which  for 
sufficient  load  may  cause  rupture.  Since  this 
experiment  was  difficult  to  carry  out,  the  same  thing 
was  at  least  partially  achieved  during  many  trials  by 
pressing  a  ball  about  a  quarter  of  an  inch  in  diameter 
onto  the  top  side  of  approximately  one  inch  square 
plates  supported  circularly  near  their  outer  edges  and 
pierced  with  a  0.01*  inch  diameter  hole.  The  behavior 
was  observed  on  a  television  screen  which  picked  up  the 
view  of  a  microscope  focussed  on  the  bottom  side. 

•  Attempt  to  use  the  Fyfe  exploding  wire  technique  at  the 
University  of  Washington  through  contracts.  Here  a 
coated  wire  in  the  bore  of  a  cylinder  of  an  inch  or 
less  in  diameter  is  exploded  by  the  controlled 
electrical  discharge  of  electricity  stored  in  a 
condenser  bank.  It  was  easy  and  very  inexpensive  to 


Oace,  Ike  0.01  nek  hole  expanded  to  3/16  nek 


before  failure,  a  atraia  of  1773%. 


28 


cause  controlled  subsurface  (to  the  hole)  rupture; 
i.e.,  complete  or  incomplete  rupture.  Thus,  an  effect 
of  anisotropy  could  be  seen.  This  was  better,  for 
purposes  of  this  report  than  plate-slap  experiments 
partly  because  rupture  was  not  preceded  by  an  enormous 
deforming  compression  of  the  material  being  tested 
before  a  tension  wave  caused  rupture.  There  was  a  very 
important  quantitative  question,  however,  about  the 
magnitude  of  the  pressure  generated,  as  well  as  a 
question  of  whether  important  high  strength  materials 
could  be  broken.  Unfortunately,  these  questions  could 
not  be  answered  before  the  work  was  terminated. 

In  the  above  cases  a  fracture  stress  is  determined  from  a 
multi-directional  stress  system  resembling  that  of  the  loaded 
crack.  Thus  the  necessity  of  a  failure  concept  for  use  of  a 
uniaxially  determined  fracture  stress  is  largely  avoided. 

In  the  case  of  the  unpressurized  crack  or  notch  deformed  in 
plane  strain,  there  is  a  stress  Y  +  YA.  acting  perpendicularly  to 
the  crack  axis  which  we  have  associated  with  the  brittle  uniaxial 
tension  test,  and  a  stress  Y L  along  the  axis.  Fracture  openings 
(preventing  crack  growth)  have  occurred  across  this  lesser 
stress,  the  material  involved  being  much  weaker  in  its  direction 
than  in  that  of  the  major  stress. 
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In  summary  this  study  relies  for  fracture  on  the  macroscopic 
stress  found  to  give  almost  nil  ductility  in  the  tension  test. 

It  is  important  to  note  that  creep  or  other  rate  effects  may 
occur  at  the  tip  of  the  crack  and,  since  the  region  to  the  point, 
x,  of  largest  stress  is  so  very  small  and  has  yielded,  that  the 
fracture  stress  is  environmentally  sensitive;  thus  progressive 
cracking  is  to  be  expected  depending  on  contaminants  and  the 
diffusion  rate  of  hydrogen  in  the  yielded  region.  The 
theoretical  fracture  stress  used  by  Griffith  is  estimated  to  be 
reached  microscopically  at  internal  flaws  but  strengths  and 
failures  at  these  individual  flaws  are  not  representative  of  the 
macroscopic  crack  strength  where  failure  covers  many  grains.  The 
situation  is  somewhat  analogous  to  that  where  it  may  be  desired 
to  determine  the  tensile  fracture  strength  of  a  structural 
component,  a  plate  pierced  by  many  holes  or  irregularities,  a 
case  where  one  measures  the  strength  by  dividing  the  breaking 
load  of  the  plate  by  its  nominal  cross-sectional  area, 
independently  of  phenomena  occurring  at  individual  holes. 

Here  as  the  load  is  increased  there  is  a  gentle  transition 
of  the  largest  stress  S,,,  from  being  proportional  to  the  load  to 
being  less  responsive  to  it  as  the  point  of  largest  stress,  x, 
moves  increasingly  sub-surface  and  finally  not  responsive  to  the 
load  at  all  when  the  stress  reaches  2.11Y  and  x  -  5p.  If  Sm  does 
reach  2.11Y  <  F,  a  marked  transition  occurs  since  from  this  point 
on  only  ductile  separation  can  occur;  the  crack  opens  up  and  the 
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yield  region  changes.  If  Sa  does  equal  F  after  the  yield  point 
is  reached  brittle  fracture  occurs.  The  failure  of  a  Charpy  Bar 
is  somewhat  similar:  above  the  transition  temperature 
corresponding  to  complete  yielding  below  the  notch  the  separation 
surfaces  open  into  a  wide  V  of  fixed  angle  which  continues  until 
final  separation  occurs;  below  the  transition  temperature  the 
separation  surfaces  remain  almost  closed  during  fracture. 

According  to  the  model  used  in  this  study,  plastic  lobes  at 
the  crack  tip  along  with  their  shear  stress  trajectories  give 
rise  to  the  maximum  stress  and  they  move  as  cracking  proceeds. 
Work  is  done  on  their  creation  and  some  of  it  is  given  back  in 
movement;  this  loss  must  be  supplied  by  the  loading.  Our  old 
attempts  to  calculate  these  lobe  generating  energies  and  losses 
were  not  completed  to  our  satisfaction;  questions  involving 
inelastic  behavior  remain.  For  a  finite  element  treatment  see 
reference  5. 

An  addition  of  hydrostatic  pressure  reduces  all  stresses  and 
thus  the  proximity  of  breakage  on  the  stress  approach  to  Kk 
failure  although  the  formulae  for  Kk,  which  may  be  derived  from 
energy  considerations,  are  independent  of  this  effect;  it  is  the 
value  for  Kk  that  differs  with  different  materials;  and  in  view 
of  the  stress  reducing  effect  of  hydrostatic  pressure,  the  value 
found  from  standard  tests  using  specimens  with  stress  free  crack 
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faces  are  not  applicable  to  pressurized  crack  faces.  We  show  how 
to  modify  the  standard  values  to  Kfcp  values  which  are  applicable. 

Finally,  let  us  examine  Griffith's  "The  Theory  of  Rupture" 
to  see  if  his  approach,  which  is  commonly  characterized  as  an 
energy  approach,  is  inconsistent  with  the  author's  approach. 

Considering  the  lack  of  appropriate  scientific  knowledge  to 
do  otherwise,  Griffith  confined  himself  to  brittle  solids  obeying 
Hooke's  Law  which  might  contain  one  or  more  cracks.  The  stresses 
in  the  solids  he  computed  from  ordinary  elasticity  theory 
including  the  stress  to  cause  rupture.  The  latter  stress  was 
taken  to  be  comparable  to  the  intrinsic  internal  pressure  of  the 
solid  (meaning  atomic  vibration  pressure  sufficient  to  cause 
vaporization)  corresponding  to  a  strain  energy  density  for 
rupture  comparable  to  the  total  heat  of  vaporization.  Rupture 
would  represent  an  enormous  temperature  rise  in  the  material  if 
the  elastic  strain  energy  did  not  go  into  surface  energy  as  in 
the  extremely  large  surfaces  associated  with  disintegration 
(e.g.,  breakage  of  very  thin  new  glass  rods)  or,  locally,  not 
concentrated  at  the  end  of  a  crack  and,  therefore,  go  into  new 
surface  if  the  crack  were  to  lengthen.  In  fact,  his  well-known 
formula  represented  a  crack  stressed  so  that  if  it  were  to 
lengthen  a  bit  (rupture)  equilibrium  would  be  maintained;  i.e., 
its  new  surface  energy  would  come  from  the  strain  energy  of  the 
solid  and  potential  energy  of  the  loads  so  that  the  total  energy 
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would  be  unchanged.  Surface  energy  density  being  a  constant, 
shorter  cracks  would  be  stable  and  longer  ones  unstable,  for  a 
fixed  loading,  since  the  stress  is  greater  for  longer  cracks. 

Thus,  we  observe  that  Griffith  relied  on  the  ideal 
theoretical  rupture  stress  for  both  cracked  and  uncracked  solids, 
but,  leaning  heavily  on  experiments  on  glass,  assumed  that  often 
cracks  were  formed  historically  which  weakened  a  solid  and  made 
its  tensile  strength  less  but  computable  for  simple  geometries 
from  the  ideal  rupture  stress  or  its  corresponding  surface  energy 
density  and  he  derived  his  failure  condition  from  an  energy 
balance. 

Thus  Griffith  would  seem  to  have  neatly  connected  both 
stress  and  energy  approaches.  But  the  crucial  thing  is  would  he 
have  used  a  different  surface  energy  density  value  than  he 
contemplated  if  he  had  thought  of  our  dodge’  of  making  the 
Griffith  derivation  trivial  by  employing  internal*1  rather  than 
with  external  loading  (and  so  avoided  the  error  in  his  first 
paper)?  We  prefer  to  think  he  would  have  done  so  and  so  agreed 
with  our  different  Kk  value  for  internal  than  external  loading. 


Making  the  pressure  cancel  the  exterior  loading  to  that  all  the  loading  ia  inside  the  crack  makes  the  derivation  of  the  Griffith  formula  extremely 
simple  since  the  energy  due  to  the  crack  is  simply  the  work  done  by  the  pressure  in  creating  volume  change  inside  the  crack  as  a  result  of  the  pressure  (see 
Appendix  0)- 
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Transverse  View  of  Mid-Thickness  Slice  Showing 


Fracture  Origin  At  Elastic-Plastic  Boundary  and  Slip  Lines 


p  =  0.005  In.  350x. 


FIGURE  1 
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TRUE  STRESS  (  PSIxlO 


from  the  data  of  (BECHTOLD)** 


True  stress  -  true  strain  curves  for 
annealed  tungsten.  Test  temp.=  523°K 


FIGURE  6 
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A  FRACTURE  MODEL:  FRACTURE  STRESS 


S 


f 


AMBIENT  PRESSURE:  p  »  0 
SET  ONE  BLOCK  ON  ANOTHER,  LOAD,  F. 
FRACTURE  (LOAD  AND)  STRESS?  Sf  «  0. 

AMBIENT  PRESSURE:  p  >  >  Y 
BLOCKS  IN  SHEATH.  LET  F/A  >  Y. 
PLASTIC  ELONGATION?  YES. 

STRESS  F/A  AT  FRACTURE?  p. 

FRACTURE  STRESS?  Sf  *  0. 


NOTE:  The  mating  surfaces  need  not  be  flat 


a)  FRACTURE  STRESS 
«  MAXIMUM  STRESS 
ACROSS  WEAKEST  PLANE, 
INCLUDING  HYDROSTATIC 
STRESS;  SAME  FOR 
NOTCHES. 


b)  LATERAL  FRACTURE 
STRESS  MAY  BE  LARGE. 


(  R.  BEEUWKES  JR.) 


19-066-1 1 59/AMC-67 


MEANING  OF  FRACTURE  STRESS 


FIGURE  7 


900  800  700  600  400 

TEMPERING  TEMPERATURE  (F) 

i  i  l  I  i 


TENSILE  STRENGTH  (  K  a  i  ) 


INFLUENCE  OF  SULFUR  ON  PLANE-STRAIN  FRACTURE  TOUGHNESS  OF  0.45C 
Ni-Cr-Mo  STEELS  (Experiments:  A.  J.  Birkle,  R.  P.  Wei,  and  G.  E. 
Pellissier,  U.  S.  Steel  Corp.  Theory:  R.  Beeuwkes,  Jr.  - 
Fracture  Stress  =2.85.5  +  Tensile  Strength;  Radii  x  10**,  in.,  ■ 

9.28,  7.83,  6.33,  4.48  for  0.008,  0.016,  0.025,  0.049  respectively.) 


FIGURE  8 
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Constant  f  Overlay 

Applicable  where  f(»F-Y)  and  p  do  not  vary 

with  Y(8»0). 

1.  Plot  K  vs  Y  data  on  single-cycle  log-log 
paper  (Keuffel  and  Esser  No.  458-100, 

No.  46-7002  or  equivalent  —  see  text)  with 
K  as  ordinate. 

2.  Position  this  overlay  so  that  its  curve 
best  fits  the  data  while  keeping  the  axes 
of  the  overlay  parallel  with  the  axes  of 
the  paper. 

3.  f  is  the  value  of  Y  at  the  point  where  the 
abscissa  *  1  line  of  the  overlay  crosses 
the  Y  axis. 

4.  f  i/wp  is  one  tenth  of  the  value  of  K 
where  the  ordinate  *  10  line  of  the 
overlay  crosses  the  K  axis. 


fc; - 

FIRST  REPORT 


FINITE  ELEMENT 
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FIGURE  1 1 


Constant  F  Overla 


Applicable  where  F  and  p  do  not  vary  with 

Y  (B=-l) . 

1.  Plot  K  vs  Y  data  on  single-cycle  log- 
log  paper  (Keuffel  and  Esser  No.  458-100, 
No.  46-7002  or  equivalent  —  see  text)  with 
K  as  ordinate. 

2.  Position  this  overlay  so  that  its  curve 
best  fits  the  data  while  keeping  the 
axes  of  the  overlay  parallel  with  the 
axes  of  the  paper. 

3.  F  is  the  value  of  Y  at  the  point  where 
the  abscissa  =  1  line  of  the  overlay 
crosses  the  Y  axis. 

4.  F  */Trp"  is  the  value  of  K  where  the  ordinate 
*  1  line  of  the  overlay  crosses  the  K  axis. 

O  F-E  Axes 
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FIGURE  12 


APPENDIX  I 


TRANSPARENT  OVERLAYS 

The  transparent  overlays  are  made  to  be  used  repeatedly 
until  worn  out  so  long  as  it  is  desired  to  determine  fracture 
stress  and  effective  crack  tip  radius  by  use  of  them  from  yield 
strength  and  corresponding  Kk  measurements. 

The  overlays  contain  a  curve,  or  curves,  each  identified  by 
a  parameter  which  conveys  their  relationship  to  a  law  relating 
fracture  stress  to  yield  strength.  The  curves  are  log-log  plots 
of  a  purely  mechanical  stress  analysis  relationship  existing 
between  KIc  yield  strength,  modulus  of  elasticity,  Poisson's 
ratio,  fracture  stress,  and  effective  crack  tip  radius,  though 
expressed  simply  in  terms  of  dimensionless  groupings  of  these 
variables  called  L  and  L  and  expressed  as  L  versus  L.  This 
relationship  is  insufficient  in  itself  to  determine  fracture 
stress  and  radius  in  the  absence  of  a  fracture  stress  law 
relating  fracture  strength  to  yield  strength  and  also  a  relation, 
here  assumed  to  be  determined  experimentally,  between  Klc  and 
yield  strength. 

In  use  the  overlay  is  laid  upon  a  log-log  plot*  of 
experimental  or  slightly  modified  experimental  K,c  versus  Y  data, 
in  such  a  way  as  to  best  match  the  data  in  all  or  part  of  its 


To  the  lame  scale  used  io  constructing  the  overlay.  Figures  114  12 


may  be  used  as  transparent  overlays  using  a  good  back  light. 
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range.  In  this  position,  the  values  of  K  and  Y  under  the  unit 
coordinate  lines  of  the  overlay  yield  the  desired  values  of 
fracture  stress  and  radius,  as  described  below.  The  match  must 
include  at  least  three  experimental  points,  for  at  least  three 
points  are  necessary  to  establish  curvature  which,  of  course, 
describes  the  shape  of  a  curve  with  a  continuous  slope. 

The  overlay  curves  themselves  are  constructed  from  a  table 
of  L  versus  A.  for  assumed  constant  values  of  the  fracture  law 
parameter  which  is  designated  here  by  0  and  covers  a  practical 
instructive  range  of  -2  £0*1  for  the  fracture  law  used  here. 

The  curves  are  graphs  of  log  l/(A.-0)  versus  log  1/[(A.-0)L]  and 
are  easily  transferred  to  transparencies  by  a  number  of  copy 
machines.*  The  curves  are  interpreted  in  terms  of  fracture 
stress  and  radius. 

Here,  as  always,  the  fracture  stress  is 


F-  Y+  YA. 


so  that 


/!  = 


F-Y 

Y 


For  example,  a  3-M  Thennofax  Copier. 
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The  fracture  stress  law  is  assumed  to  be  linear  in  the  Kk,  Y 
region  being  matched;  i.e., 

F-  (1  +  P)  Y+f 

so  that 

v  Y  Y  Y 

and ,  hence ,  A.  -  P  =  f/Y 

The  abscissa  on  the  transparency  is 

log  1/  (A.-P)  =log  Y/  f -log  Y-logf 

so  that  Y  =  f  where  Y/f  =  1  and  log  Y/f  =  log  1 /(A  -  0)  =0. 

The  ordinate  is 

log  1/KA.-P)  L]  =  log  [(Y/f)  (K/Yjnp)] 

=  logK’-logfVrcp 

and 

K  =  fyj tip  where  1/  [  (A.-P)  L]  =1 
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Thus,  if  one  of  the  /S  curves  of  the  transparent  overlay  can 
be  matched  to  a  log-log  plot  of  experimental  values  of  Kfc  versus 
Y,  the  value  of  the  Y  intersected  by  the  unit  coordinate  line  of 
the  overlay  will  be  Y  =  f,  for  what  is  log  Y  on  the  plot  is  [log 
Y  -  log  f];  i.e.,  0  on  the  overlay  when  it  is  in  the  matching 
position.  Correspondingly  and  similarly,  the  horizontal  unit 

coordinate  line  of  the  overlay  is  where  K  =  fjnp  and  from  this 

v/np  is  derived  by  dividing  out  the  value  of  f.  Therefore,  yfrt p 
as  well  as  the  constants  /S  and  f  of  our  fracture  law  are  known. 

The  overlays  and  graphs  for  the  examples  in  the  First  Report 
were  for  0  =  0;  i.e.,  F  =  Y  +  f,  and  j8  =  -1  (i.e.,  F  =  f, 
independent  of  Y) .  The  curves  for  making  these  overlays  and  the 
experimental  data  were  plotted  on  the  largest  size  single  cycle 
log- log  paper  locally  available  that  would  fit  into  this  report 
with  K  on  the  vertical  and  Y  on  the  horizontal  scale.  This  size 
paper  and  the  corresponding  overlays  are  considered  so  generally 
useful  that  such  overlays  are  included  in  this  report  and  are 
called  the  f  (for  /3  =  0)  and  F  (for  /3  =  -1)  overlays. 

However,  to  make  an  overlay  to  cover  the  larger  range  of 
experimental  variables  and  ranges  of  /S  that  may  be  encountered,  2 
x  3  log-log  paper,  whose  size  also  conforms  to  that  of  the  report 
had  to  be  used.  Unlike  the  scale  of  the  single  cyle  paper  above, 
this  selection  of  paper  means  that  for  this  wider  range  the  scale 
is  such  that  the  Y  values  (plotted  on  the  horizontal  axis,  as 
before)  increase  from  right  to  left. 
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APPENDIX  II 


ENERGY  RELATIONS  IN  SUPERPOSITION  OP  HYDROSTATIC  PRES8URE 

In  view  of  the  fact  that,  in  contrast  to  three-dimensional 
problems,  in  plane  strain  a  lateral  stress  is  required  to  keep  the 
thickness  constant,  and  in  plane  stress  there  is  no  lateral  stress, 
it  may  be  wondered  just  how  the  superposition  of  hydrostatic 
pressure  is  to  be  carried  out  (whether  lateral  pressure  is  to  be 
employed)  .  Here  the  superposition  is  taken  in  such  a  way  as  to 
preserve  the  plane  strain  and  stress  conditions.  Thus,  to  modify 
the  general  expression  for  energy  per  unit  volume  in  terms  of  the 
stresses  for  the  inclusion  of  hydrostatic  pressure,  p: 

•  for  three  dimensions  add  -  p  to  each  stress 

•  for  plane  stress  make  the  stress  in  the  thickness 

direction  zero  and  add  nothing  to  it 

•  for  plane  strain  make  the  stress  in  the  thickness 

direction  a  (ax  +  ay)  in  which  a  without  a  subscript  is 
Poissons 's  ratio,  and  add  p  to  ax  and  to  ay,  including  the 
ax  and  ay  in  a(ax  +  ay)  . 

Each  modified  expression  for  energy  is  the  integrand  of  an 
integral  over  the  volume  of  a  body  for  the  total  energy  of  that 
body;  the  part,  in  each  case,  without  p  is  the  whole  energy  before 
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superposition  of  p  and  the  part  containing  terms  in  p  is  the 
decrease  in  energy  caused  by  superposition  of  pressure  p.  The 
latter  consists  of  a  simple  term  in  p2  and  a  part  in  which  the  sum 
of  the  original  stresses  is  multiplied  by  p.  By  replacing  each 
such  stress  by  its  isotropic  Hooke's  Law  expression  of  the  type  a 
=  XA  +  2 he,  it  is  seen  that  the  sum  represents  the  unit  dilatation 
so  that  the  integral  of  this  part  is  a  constant  multiplied  by 
pfAdV,  i.e.,  pressure  multiplied  by  the  whole  dilatation  caused  by 
the  original  stresses. 

This  dilatation  was  produced  by  an  original  exterior  boundary 
loading  stress  +  p  since  this  is  the  p  that  is  being  eliminated  on 
the  exterior  boundary  by  the  hydrostatic  addition  of  -p.  The  total 
dilatation  jAdV  above  is  expressed  in  terms  of  the  original  loading 
stress  +p  by  means  of  a  theorem  in  Love's  "The  Mathematical  Theory 
of  Elasticity",  Article  123. 

It  is  found  that  (since  the  crack  has  no  volume  in  the 
unstrained  condition)  the  total  energy  contribution  caused  by  the 
hydrostatic  -p  addition  is: 

•  for  three  dimensions  -p2V/2K 

Where  the  bulk  modulus  K  =  E/[3(l  -  2a)] 

•  for  plane  stress  -p2  V/2K, 

where  the  bulk  area  modulus  K,  =  E/[2(l  -  a)] 

•  for  plane  strain  -p2V/2K2 

Where  the  bulk  area  modulus  K2  =  E/[2(l  -  a)  (1  -  2a)) 
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as  though  there  was  no  crack  at  all. 


Where  hydrostatic  pressure  cancels  the  external  load  leaving 
the  crack  under  internal  pressure,  these  simple  terms  are  the 
amounts  the  total  original  energy  has  decreased;  what  energy 
remains,  since  the  only  load  is  now  in  the  crack,  is  due  to  that 
pressure  in  the  crack.  This  energy,  because  displacement  in  the 
crack  is  proportional  to  the  internal  pressure,  is  the  pressure 
multiplied  by  half  the  volume  it  produces. 

Now,  consider  from  this  background  of  superimposed  hydrostatic 
stress,  Griffith's  analytical  difficulty.  He  states,  "A  solution 
of  this  problem  was  given  in  a  paper  read  in  1920  but  in  the 
solution  there  given  the  calculation  of  the  strain  ern '  gy  was 
erroneous  in  that  the  expression  used  for  the  stresses  gave  values 
at  infinity  differing  from  the  postulated  uniform  stress  at 
infinity  by  an  amount  which,  though  infinitesimal,  yet  made  a 
finite  contribution  to  the  energy  when  intregated  around  the 
infinite  boundary.  This  difficulty  has  been  overcome  by  slightly 
modifying  the  expressions  for  the  stresses,  so  as  to  make  this 
contribution  to  the  energy  vanish,  and  the  corrected  condition  for 
the  rupture  of  a  thin  cracked  plate  under  a  uniform  pull  applied  in 

its  plane  at  its  outer  edge  is  R~^2Et/nc  where...".4 

Now  suppose  we  have  a  plate  containing  a  crack  (or  a  hole  for 
that  matter)  under  internal  pressure  p  only,  superposition  of  a 
hydrostatic  tension  p  on  that  plate  will  eliminate  the  internal 
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pressure  p  and  load  the  plate  on  its  outer  edge  with  a  tension 
stress  acting  normally  to  its  outer  edge;  there  are  no  secondary 
stresses  even  though  no  limit  has  been  put  on  the  size  or  shape  of 
the  plate.  Any  energy  input  caused  by  the  presence  of  a  crack  is 
due  to  the  action  of  the  original  internal  pressure  on  the  sides  of 
the  crack.'  Within  reason,  it  would  seem  that  little  error  in  this 
energy  input  would  be  made  if  the  maximum  opening  of  the  crack  is 
considered  known,  if  the  crack  is  assumed  to  have  gone  into  an 
elliptical  shape  even  where  this  is  not  the  case;  e.g.,  if  the 
cracked  plate  is  short  or  narrow  in  comparison  to  the  size  of  the 
crack. 


54 


REFERENCES 


1.  BEEUWKES,  R.,  Jr.  "Analysis  of  Failures,"  Proceedings  of  the  Third 
Sagamore  Ordnance  Materials  Research  Conference,  Analysis  of  Failures, 

Duke  University,  5-7  December  1956,  p.  89-163.  Printed  by  Syracuse  Uni¬ 
versity  Research  Institute  (Correction:  Multiply  (N  -  Nq)  by  (3/4)*/3 

2.  BEEUWKES,  R.,  Jr.  "Characteristics  of  Crack  Failure,"  Surfaces  and 
Interfaces,  v.  II,  Syracuse  University  Press,  Syracuse,  NY,  1968. 

3.  BEEUWKES,  R.,  Jr.  "On  Status  of  Correlation  of  Crack  Failure  with  Mate¬ 
rial  Properties  in  Metals,"  U.S.  Army  Research  Laboratory,  AMMRC  TN  69-01, 
March  1969. 

4.  BEEUWKES,  R. ,  Jr.,  HANSEN,  A.  M.,  and  TREMBLAY,  R.  J.  "Determination 

of  Fracture  Stress  and  Effective  Crack  Tip  Radius  from  Toughness  (Kjc)  an<* 
Yield  Strength  (Y),"  U.S.  Army  Research  Laboratory,  AMMRC  TR  78-44, 

October  1978. 

5.  KOBAYASHI,  A.  S.,  CHIU,  S.  T.,  and  BEEUWKES,  R. ,  Jr.  "Elastic-Plastic 
State  in  a  Plate  With  an  Extending  Crack,"  Proceedings  of  the  Army 
Symposium  on  Solid  Mechanics  -  1970,  Lightweight  Structures:  also  U.S. 

Army  Research  Laboratory,  AMMRC  MS  70-5,  v.  IV,  December  1970. 

6.  KOBAYASHI,  A.  S.,  CHIU,  S.  T. ,  and  BEEUWKES,  R. ,  Jr.  "Elastic-Plastic 
State  and  Fracture  of  a  Blunt  Notched  Tension  Plate,"  in  Mechanical 
Behavior  of  Materials.  The  Society  of  Materials  Science,  v.  1,  Japan, 

1972,  p.  446-457. 

7.  TRACEY,  D.  M. ,  and  FREESE,  C.  E.  "Plasticity  Near  a  Blunt  Flaw  Under  Remote 
Tension,"  ASTM  STP  995,  1989,  p.  93-106. 

8.  TRACEY,  D.  M. ,  and  FREESE,  C.  E.  "Importance  of  Crack  Tip  Shape  in 
Elastic-Plastic  Fracture  Analysis,"  Proceedings  of  the  Second  Army  Confer¬ 
ence  on  Applied  Math  and  Computing,"  1985,  p.  359-371,  U.S.  Army  Materials 
and  Mechanics  Research  Center,  Watertown,  MA. 

9.  KE,  J.  S. ,  and  LIU,  H.  W.  "Thickness  Effect  on  Crack  Tip  Deformation  at 
Fracture,"  Syracuse  University,  Contract  Number  DAAG46-73-C-0209 ,  Final 
Report,  AMMRC  CTR  74-57,  Volume  IV,  December  1976. 

10.  KE,  J.  S.,  and  LIU,  H.  W.  "The  Measurements  of  Fracture  Toughness  of 
Ductile  Materials,"  Presented  at  the  Symposium  on  Fracture  and  Fatigue, 
George  Washington  University,  Washington,  DC.  Pergamon  Press,  Oxford,  NJ, 
1972. 

11.  GRIFFITH,  A.  A.  "The  Theory  of  Rupture,"  Proceedings  of  the  International 
Congress  for  Applied  Mechanics,  Delft,  Holland,  1924. 

12.  GRIFFITH,  A.  A.  "The  Phenomena  of  Rupture  and  Flow  in  Solids,"  Phil. 

Trans.  Roy.  Soc .  London,  A.,  v.  221,  1921,  p.  163-198. 

13.  TAGGART,  R. ,  and  WAHI,  K.  K.  "Crack  Opening  Displacement  During  the 
Tensile  Loading  of  Ductile  and  Brittle  Notched  Plates,"  University  of 
Washington,  Seattle,  under  subcontract  to  Syracuse  University  Prime 
Contract  DA-31-124-ARO(D)-l 12 ,  Final  Report,  AMMRC  CTR  72-33,  February 
1972;  also,  abridged  version  entitled  "The  Relationship  between  the 
Fracture  Toughness  and  the  Crack  Tip  Radius,"  by  TAGGART,  R. ,  WAHI,  K.  K. , 
and  BEEUWKES,  R.,  Jr.,  ASTM  STP  605,  1976,  p.  62. 


55 


14.  GARDE,  A.,  and  WEISS,  V.  "Brittle  Fracture  Initiation  at  the  Elastic- 
Plastic  Interface,"  Syracuse  University  Research  Institute,  Contract 
DA-31 -124-ARO(D)-l 12 ,  Interim  Technical  Report,  AMMRC-CR  70-29, 

February  1970;  also,  "Brittle  Crack  Initiation  at  the  Elastic-Plastic 
Interface,"  Met.  Trans.  AIME,  v.  3,  November  1972,  p.  2811. 

15.  BEEUWKES,  R. ,  Jr.  "Stresses  in  Thick  Walled  Cylinders,"  presented  at  the 
International  Congress  for  Applied  Mechanics,  Paris,  France,  10  October 
1946.  Watertown  Arsenal  Laboratory  Report  WAL  730/419,  1946. 

16.  SMITH,  J.  H.  "Some  Experiments  on  Fatigue  of  Metals,"  Journal  of  the 
Iron  and  Steel  Institute,  September  1910,  p.  246-318. 

17.  SMITH,  J.  H.,  and  WEDGEWOOD,  G.  A.  "Stress-Strain  Loops  for  Steel  in 
the  Cyclic  State,"  Journal  of  the  Iron  and  Steel  Institute,  1915, 

p.  365-397. 

18.  BATSON,  R.  G.,  and  HYDE,  J.  H.  "Mechanical  Testing,"  The  D.  U.  Technical 
Series,  published  by  E.  P.  Dutton,  p.  224-226 

19.  BECHTOLD,  J.  H. ,  and  SHEWMON,  P.  G.  "Flow  and  Fracture  Characteristics 
of  Annealed  Tungsten,"  Trans  ASM,  v.  46,  1954,  p.  403  (See  Figure  8). 


No.  of 
Copies 


DISTRIBUTION  LIST 


To 


Director,  U.S.  Army  Research  Laboratory,  2800  Powder  Mill  Road,  Adelphi, 

MD  20783-1197 

1  ATTN:  AMSRL-OP-CI-AD,  Technical  Publishing  Branch 

1  AMSRL-OP-CI-AD,  Records  Management  Administrator 

Commander,  Defense  Technical  Information  Center,  Cameron  Station,  Building  5, 
5010  Duke  Street,  Alexandria,  VA  23304-6145 

2  ATTN:  DTIC-FDAC 

Commander,  U.S.  Army  Materiel  Command,  5001  Eisenhower  Avenue,  Alexandna, 
VA  22333 

1  ATTN:  AMCSCI 

Director,  U.S.  Army  Research  Laboratory,  Watertown,  MA  02172-0001 

2  ATTN:  AMSRL-OP-WT-IS,  Technical  Library 

5  Author 


57 


